Universal features of the defect-induced damping of lattice vibrations 



A. Cano,i'EI A. P. Lcvanyuk,2.[] and S. A. Minyukov^'H 

^ Instituto de Astrofisica de Andalucia CSIC, PO Box 3004, E-18080 Granada, Spain 
^ Departamento de Fisica de la Materia Condensada C-III, Universidad Autonoma de Madrid, E-28049 Madrid, Spain 
'^Institute of Crystallography, Russian Academy of Sciences, Leninskii Prospect 59, Moscow 117333, Russia 

(Dated: February 5, 2008) 

It is shown that any defect gives an Ohmic contribution to the damping of any normal mode of the 
crystal lattice with nonzero wavevector which does not vanish at zero temperature. This explains 
the large phason damping observed at low temperatures in incommensurate phases, and might be 
a key factor to understand the linear-in-T specific heat observed in a number of real dielectrics at 
low enough temperatures. 



It is well known since long ago that the defect-induced 
damping of the normal modes of the crystal lattice is a 
key factor to understand the real dielectric losses , at- 
tenuation of sound 0, spin lattice relaxation rates, etc. 
This is especially clear at low temperatures, where the 
real damping cannot have its origin in the well studied 
effects of anharmonism: a strong dependence on temper- 
ature ~ -f- is predicted for the phonon attenuation 
01 which is incompatible with the low-tcmpcrature 
data 0. Recently it has been argued that, at low tem- 
peratures, the defect-induced damping can play an even 
more far reaching role as long as it may be the reason 
of the deviations from the Debye law observed in the 
low-temperature specific heat of real dielectrics 0, 0] . In 
this paper we show that the defect-induced damping is 
unexpectedly universal in the sense that any defect gives 
an Ohmic contribution to the damping of any normal 
mode of the crystal lattice with nonzero wavevector. This 
feature reinforces the candidature of the defect-induced 
damping as the main reason for the aforementioned de- 
viations from the Debye law. 

In accordance with previous results restricted to a 
k = polar mode, the presence of a random distribu- 
tion of charge defects makes the damping of this mode 
to be Ohmic at low frequencies, i.e., without frequency 
dispersion 0. This can be easily understood as follows. 
When applying a homogeneous electric field, i.e., the 
force conjugated to the homogeneous polarization, the 
charge defects do displace. These displacements are ac- 
companied with (local) displacements of the correspond- 
ing lattice. So, when the field varies in time, the motion 
of the charges provokes the radiation of acoustic waves 
whose intensity is cx lu'^ (not cx as it could be inferred 
from Rcf. for the case of local deformations). This 
radiation gives rise to dielectric losses. In Ref. |3 it was 
shown that, in the case of a correlated density of charge 
defects, this result is also valid for most of the modes in 
a polar branch. In the following we shall show that, for 
any lattice mode with nonzero wavevector, any type of 
defect will also induce the same type of damping. 

For the sake of illustration, let us consider a polar 
branch as in the previous papers. It is worth mentioning. 



first of all, that the presence of defects means that there 
is a redistribution of charges with respect to that in the 
ideal lattice. This allows us to model any defect as a set 
of electric multipoles. As we shall see, an inhomogeneous 
electric field varying in time has the above effect on any 
electric multipole. That is, the field induces the radia- 
tion of acoustic waves through the multipoles and, conse- 
quently, gives rise to losses. To be specific, let us consider 
randomly distributed dipolar defects. The charge density 
associated with these defects can be written as 



p(r) = ^(p,.V)<5(r-r,), 



(1) 



where p, = qdi and arc (random) vectors representing 
the dipolar momenta of the defects and their positions, 
respectively. (For the sake of simplicity, we further con- 
sider that the dipoles are oriented randomly.) 

The frequency dependence of the damping constant 7 
of a polar mode can be obtained, in the low-frequency 
regime we are interested in, from the frequency depen- 
dence of the imaginary part of the corresponding suscep- 
tibility x(w,k): 7(lj ^ 0,k) cx uj-'^x"{^ ^ 0,k). To 
compute such a susceptibility, we first take into account 
that, in Fourier space, the polarization can be written as 

P(w,k) =xo(w,k)E(c^,k) 
1 



V 



{p(r)u{r,t))e-'^^-''-'^''>drdt. (2) 



Here xo is the (scalar) susceptibility in absence of defects 
(the medium is assumed to be isotropic), E is the (inho- 
mogeneous) electric field, V is the volume of the system, 
p is the charge distribution due to the defects, u is the 
acoustic displacement vector, and (...) denotes statisti- 
cal average. The equation of motion for the displacement 
vector is 



ii = Ct V^u + (cf - )V(V • u) + pE/^, 



(3) 



where c; and ct are the longitudinal and the transversal 
velocities of sound, respectively, and p, the density of the 
medium. [The last term in this equation is just the force 
that the electric field exerts on the lattice due to the 
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presence of the (charges) defects.] In Fourier space, the 
i-th component of this displacement vector then can be 
written as 



u,;(w,k) ^ -A./(tj,k) 



V 



^p(k-k')i?,(c.,k'), 



(4) 



where 



Substituting Eq. Q into Eq. (0) we find that 



V r dk' 



A-(k-k')A,r.i(c.,k')£;,(c^,k), 

(6) 



where K{\^) is defined from the "spectral density" of the 
charge distribution: (p(k)p(k')) ^ K{\^)V-^5{\i + 'k') 0. 
In accordance with Eq. (Q, this magnitude is 



(7) 



in our case, where N is the density of defects. 

As long as we are considering an isotropic medium, 
after integrating over wavevectors we find that the con- 
tribution of the non-diagonal terms of A^^^ in Eq. ® is 
zero. This means that the contribution to the suscepti- 
bility due to the defects can be written as 



Xdcf(t^,k) 



V 



M J (27r)3 



if(k-k')A,:-i(c.,k') (8) 



(here double subscript docs not imply summation). The 
main contribution to the integral in Eq. JHJ comes from 
the poles of k-^{Lo,\i!). We arc interested in the case 
of small frequencies. Taking into account that the poles 
of A,^"'^((jj, k') are such that k' ^ oj/c {ci,ct ^ c), in the 
integrand of Eq. © we can take 

K{k-k') ^ Np'^k^/{3V), (9) 

for large wavevectors (fc i^/c). We then have 

// / „ , , N N p^k^ uj , , 

iVat UJD 

where A^at = ^at'^ is the atomic density and ujd = c/lat 
is the Debye frequency. In this expression we see that, 
though the damping vanishes for wavevectors tending to 
zero, it is quite significant for wavevectors close to the 
boundary of the corresponding Brillouin zone (where the 
factor pk is ~ g) . 

At this point, it is worth mentioning that, at large 
enough wavevectors, the modes of any branch are effec- 
tively polar (strictly speaking, the label polar or nonpolar 



refers to k = 0). Therefore, the above considered mech- 
anism for damping gives a contribution for either optic 
or acoustic modes. In the latter case, the polarization 
may have its origin in either ions or electrons. It is worth 
noticing that, in any case, the strength of the damping is 
the same for a; — > [the only parameters entering in Eq. 
are elastic constants and the multipole moment]. 
It is also worth mentioning that similar results are ob- 
tained not only for higher order electric multipoles, what 
can be easily seen by following similar arguments, but 
also in the case of multipoles of different nature. In 
fact, the change in the charge distribution of a crystal 
is only one of the several aspects of the infiuence of de- 
fects. Among these other aspects, much attention has 
been paid to the elastic one for which, the modelling of 
the defects as "force multipoles," has been proved to be 
quite useful (see, e.g., Ref. 0]). Inhomogeneous strains 
and stresses then play a role analogous to that of the in- 
homogeneous polarization and electric fields considered 
above. 

Defects can also be considered as more general mul- 
tipoles, say "optic multipoles." This became apparent 
when studying phase transitions where, acting as a (lo- 
cal) field conjugated with the corresponding order pa- 
rameter, defects may play a key role (see, e.g., Ref. |ll|). 
In view of this, the above mechanisms may explain, for 
instance, the unexpected phason damping observed at 
low temperatures (see, e.g., Ref. 0]). The theory of an- 
harmonic crystals cannot explain this damping as long as 
this theory predicts a very strong dependence on temper- 
ature 0]. But taking into account that, in terms of the 
normal phase, phasons are optic or acoustic vibrations 
with large wavevectors (close to that of the correspond- 
ing incommensurate structure), it is reasonable to expect 
a significant contribution due to the above mechanism of 
damping. 

As an interesting consequence of the above results, it 
is the fact that a damping with such features provides 
a natural explanation for the, otherwise difficult to ex- 
plain, (quasi)linear-in-r contributions observed in the 
low-temperature specific heat of a number of incommen- 
surate dielectrics and charge-density- wave systems [l^j . 
The mechanism of damping we have discussed naturally 
explains also the non-Debye specific heat observed in a 
number of ferroelectrics some time ago 0], which was 
cxperimctally linked to the presence of defects Let 
us detail this point to some extent. 

In principle, the low-temperature properties of di- 
electrics are determined by the corresponding density of 
(pho non) states at frequencies tending to zero (see, e.g., 
Ref. [13 )■ In real dielectrics, i.e., dielectrics containing 
defects, the general features of such a density of states, 
however, are normally unknown. Whether this density 
of states remains finite at a; = is a question of particu- 
lar interest: In this case, the crystal will exhibit, among 
other characteristic features, a lincar-in-T specific heat at 
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low enough temperatures; what is normally considered as 
a fingerprint of glasses (see, e.g., Ref. 0]). 

When the concentration of defects is small, it seems 
plausible to study the low-temperature properties of real 
systems, instead of focusing on the corresponding den- 
sity of states, by following the approach exposed in Refs. 

(see also Ref. 0). When doing so, one first consid- 
ers the phononic scenario of the clean (ideal) crystal and 
then encapsulates the influence of defects into a defect- 
induced damping of the corresponding phonons 0, 0|. 
The possibility of identifying a reservoir where the energy 
lost by the damped phonons can be transferred to, with- 
out altering the properties of such a reservoir, allows one 
to further compute the corresponding low-temperature 
properties within the already developed framework of 
quantum dissipative systems (see, e.g., Ref. ^3)- Long- 
wavelength acoustic phonons brings about the mentioned 
reservoir , and the validity of the approach is guaran- 
teed by virtue of the smallness of the defect concentration 
— i.e., the weakness of coupling between the (sub)system 
of interest and the reservoir. Within this context, to 
have an Ohmic damping for an harmonic oscillator is 
equivalent to have a finite density of states at w = 
and, though it has not been demonstrated explicitly, it is 
quite expectable that the same happens when one con- 
siders a whole crystal (i.e., not a single oscillator but a 
set of them). Bearing in mind that the above Ohmic 
damping is precisely the one we have demonstrated to be 
"universal," the specific heat of real dielectrics has to be 
expected to be linear-in-T at low enough temperatures. 

A natural question is the relevance of these results for 
glasses. Strictly speaking, glasses are beyond of the ap- 
plicability of the approach in Refs. as long as this is 
designed to deal with small concentrations of defects only. 
However, it is quite reasonable to speculate that, behind 
the low-temperature specific features of glasses, there are 
significant contributions that may be interpreted as due 
to the damping of the corresponding vibrational excita- 
tions. 

In summary, we have shown that any defect gives an 
Ohmic contribution to the damping of any lattice vibra- 
tion with nonzero wavevector. This contribution does not 
vanish at zero temperature and, consequently, may be a 
key point to understand the low-temperature properties 
of real systems that, despite a small concentration of de- 
fects, exhibit universal glass-like features. 
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